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Abstract. We review the problem of proton beam bombardment of solar chromosphere 
considering the self-consistent NLTE polarized radiation transfer in hydrogen lines. Several 
observations indicate a linear polarization of the Ha line of the order of 5% or higher and 
preferentially in radial direction. This polarization is often explained as anisotropic colli- 
sional excitation of the n = 3 level by vertical proton beams. Our calculations indicate that 
deceleration of the proton beam with initial power-law energy distribution together with 
increased electron and proton densities in the Ha forming layers lead to a negligible line 
polarization. Thus the proton beams seem not to be a good candidate for explanation of the 
observed polarization degree. On the other hand, the effect of electric return currents could 
perhaps provide a better explanation of the observed linear polarization. We report the new 
calculations of this effect. 
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1. Introduction 

There is an observational evidence for the fast 
electron beams (10-100 keV) and also for the 
fast proton beams (~ 10 MeV) bombarding 
the solar chromosphere during the i mpulsive 
and gradual phase of solar flares dKorchakl 
Il967t lOrrall & Zirkerl [l976h . This bombard- 
ment is consistent with the "standard model" 
of solar flare which assumes an injection of 
the high energetic particles from a coronal 
reconnection site to the chromosphere. It re- 
sults in heating and nonthermal excitation of 
the chromospheric gas. The standard semiem- 
pirical models of th e flaring chromosphere 
(Machado et al ] ll980h are based on the series 



of continuum and lines observations and do not 
take into account a nonthermal excitation; they 
rather overestimates the chromospheric tem- 
perature to explain an increased radiative emis- 
sion. 

Also the low energy proton beams (below 
1 MeV) could play a significant role in the 
flare physics but their presence in the chro- 
mospheric layers is still uncertain due to their 
negligible bremsstrahlung radiation. There are 
however different ways how to diagnose them. 
For instance: 

1. Lya red wing emission due to charge 
exchange effect dCanfield & Changl 119851: 
IZhaoetalJll998l) . 
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2. Modification of the lines intensity 
profiles due to nonthermal exc itation 
ilHenoux e7aill993tlXu et alj|2005ah . 

3. Detection of a linear polarization of the 
lines due to an anisotropic collisional ex- 
citation (i.e. impact atomic polarization). 

Moreover, several unanswered questions still 
remain about the electron beams. Especially 
the physics of the return currents creation and 
their influence on lines formation. 

Measurement of linear polarization in so- 
lar flares is a complicated task due to the high 
spatial and time gradients. Lot of observations 
report the polarizati on degree of Ha of the or- 
der of few percent dHenoux & Chambe 1990t 
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Vogt & Henouxll999HHanaokal2003tlXu et al.1 
2005b). Direction of this polarization is usually 
found to be either radial or tangential with re- 
spect to the limb and its degree is up to 5% 
or even exceeds 10% in the few cases. The 
standard interpretation of these measurements 
is anisotropic excitation by proton and/or elec- 
tron beams. By contrast, no measurable linear 
polarization above the noise level at 0.1% was 
found in the s et of many different flares by an- 
other authors dBianda et al.ll2003l [2005b . 

Our objective is to find out if the measure- 
ments of a linear polarization degree of the hy- 
drogen Ha line are a promising tool for the 
diagnostics of the low-energy proton beams. 
Moreover, we present a first quantitative es- 
timation of the polarizing effect of the return 
electric currents taking properly into account 
depolarizing collisions and radiative transfer. 

We start with the standard semiempir- 
ical model of the flaring atmosphere Fl 
dMachado et al.1 1 19801) with a fixed tempera- 
ture structure. We include the nonthermal col- 
lisional rates of excitation and ionization to ob- 
tain the new line profiles. We use a stationary 
modeling of the polarized radiation transfer in 
the chromospheric hydrogen lines using a mul- 
tilevel NLTE polarized transfer code. 

2. Ha impact polarization 

Let us assume a beam of the unidirectional 
charged colliders and let their energy be high 
enough to excite the n — 3 level (i.e. higher 
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Fig. 1. Linear polarization of the photons emitted in 
the Ho- line perpendicularly to the collisional trajec- 
tory (i.e. in the direction of a maximum polarization 
degree). For both the electrons and protons impacts, 
low energies of collisions lead to linear polarization 
of the emitted photons in the direction parallel to the 
collisional trajectory. A polarization direction above 
the so-called turnover energy is perpendicular to it. 



than 12.1 eV). In fact, the n — 3 level is com- 
posed of five quasi-degenerated fine structure 
levels. Some of them can be polarized, hence 
the photons emitted in the Ha (and Ly/J) line 
are in general polarized. An anisotropical exci- 
tation of these levels by either protons or elec- 
trons is one of the processes which can lead 
to polarization of these levels. The orientation 
of electric vector of this radiation can be either 
parallel or perpendicular to the velocity of the 
beam (cf. FigureQ]). A beam is expected to fol- 
low the direction of magnetic field which is as- 
sumed to be vertical in the atmosphere; hence 
the linear polarization is either radial or tan- 
gential with respect to the solar limb. 

The assumptions above should lead to an 
upper theoretical limit for the Ha linear po- 
larization. The coordinate system definition is 
given in Figure [2] In such a reference frame, 
only the Stokes parameters / and Q are in gen- 
eral nonzero. We are interested especially in a 
modeling of the radial polarization. It follows 
from Figure Q] that this orientation is associ- 
ated with the low-energy electrons or the low- 
energy protons moving vertically in the atmo- 
sphere 



1 The case of horizontal motions at high energies 
is unlikely for the physical reasons. 



2 



Stepan: Polarization diagnostics of proton beams in solar flares 



Particle trajectory 



*/ o 



X 



Fig. 2. We use a right-handed coordinate system 
with /-axis in the collisional plane and Z-axis in the 
line of sight direction. Due to cylindrical symmetry 
of the model, only Stokes parameters / and Q do not 
generally vanish in this coordinate system. 



3. Formalism and methods 

A vertical magnetic field is assum ed to be of 
the st rength of few hundreds gauss dVogt et al.1 
Il997l) . It can be shown that in this regime and 
in the plasma conditions under consideration 
(see below) fine structure splitting of the hy- 
drogen levels up to n — 3 is a go od approxi- 
mation (Sahal- Brechot et ail l996). The hyper- 
fine splitting of the levels can be neglected be- 
cause it does not s ignificantly affect a lin ear 
polarization degree dBommier et a l] |l986al) . In 
addition, a lifetime of the fine structure levels 
is reduced by collisions with the background 
protons and electrons and thus the hyperfine 
levels of the excited states completely over- 
lap. To describe an atomic state we adopt the 
formalism of the atomic density matrix in the 
basis of irreducible tensorial operators p^(nl j) 
(lFanolll957h . with the traditional meaning of 
the symbols. We neglect all the quantum co- 
herences between different fine structure levels 
due to their large separation in comparison to 
their width and due to the selection rules for 
the optical transitions. A magnetic field is able 
to destroy all the quantum coherences between 
the Zeeman sublevels of any level. In the for- 



malism of irreducible tensors it means that all 
the multipole components of p with Q + 
are identically zero, pg(nlj) = 6Qopg(nlj). 
On the other hand, the strength of a magnetic 
field is not so high to induce a Zeeman split- 
ting large enough to lead to the complicated 
level-crossing effects. Finally, because of the 
assumption of absence of any circularly polar- 
ized radiation, the odd ranks K of the density 
matrix are identically zero. All these assump- 
tions lead to the simplifications of the formal- 
ism. 

The results of our modeling are strongly 
dependent on the collisional cross-sections of 
the different excitation and ionization transi- 
tions. We consider the following collisional 
processes and cross-sections data: 



- Excitation and charge exchange for a 
proton beam using the close-coupling 
cross- sections data of iBalanca & Feautrierl 
(119981) . 

- The dipolar transitions between the fine 
structure levels nlj — > nl + 1 / for 
interaction with the background elec- 
trons, protons, and the beam are calcu- 
lated using the semiclas sical theory of 
ISahal-Brechot et all dl996l) . 

- Excitation of the levels populations by the 
background electrons (the cross-sections 
data are taken from the AMDIS database, 
http://www-amdis.iaea.org). 



For the purposes of impact polarization stud- 
ies it is necessary to calculate all the collisional 
transition rates 



C 



nlj^n'l'j' 



j d 3 v/(v)v^, ; ,.,(v) 



(1) 



for transitions p^(nlj) — > pF («'/'/). These 
rates enter the equations of statistical equi- 
librium together with the radiative rates 
R K i^ K ',,,; in the so-called impact ap- 
proximation ( Landi Degl' Innocent! Il984t 
iBommier & Sahal-Brechoj Il99ll) . The equa- 
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tions of statistical equilibrium have the 
forrrQ 

+ R n»Zn'M^ = • ® 

nljK 

To take properly into account a coupling of 
the atomic states at the different atmospheric 
points, it is necessary to solve a system of 
equations © coupled by the radiative transfer 
equation for Stokes vector / = (/, Q) T . Its evo- 
lution along the path s is given by 

^-=J-K1. (3) 

as 

In this equation, J is the emissivity vector and 
K is the so-called propagat ion m atrix. We use 
our multigrid code (Stepan 2006) to solve this 
problem in the plan-parallel geometry, out of 
the LTE approximation (the case of the chro- 
mospheric Ha line). 

4. Results for the proton beams 

The temperature structure of the atmosphere is 
fixed and given by the Fl model. A new ion- 
ization degree and the atomic states along the 
chromosphere are obtained by solution of the 
system of the Eqs. (O and The initial en- 
ergy distribution of the proton beam at the top 
of the chromosphere is (as usually) expected to 
be given by the power-law 

F(E) ~ E~ s , E> E , 

with the lower energy cut-off E c and the spec- 
tral index 6. After crossing some column depth, 
the distribution is modified by collisions with 
various atmospheric species (see Figu re [3]), 
(lEmsiid 119781: ICanfield & Chang! 119851) . The 
protons most efficient in creation of the n — 3 
level p olarization have a small e nergy about 
5 keV dBalanca & Feautrienfl998h . A number 
of such protons in the region of the Ha line 
core formation is very small compared to the 
number of proton s of an energy about E c a t the 
injection site (cf. IVogt et al.lll997l 120011) . As 

2 Additional terms have to be added to take into 
account the bound-free transitions. They are not ex- 
pressed in the Eq. 10 for the reasons of simplicity. 
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Fig. 3. Energy distribution of a beam flux (in arbi- 
trary units), S = 4, £ c = 150 keV. The initial sharp- 
peaked power-law energy distribution (upper curve) 
becomes flatter in the region where the slowest pro- 
tons are stopped (lower curve). 
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Fig. 4. Electron density for the model Fl. The plot 
corresponds to the case of 6 = 4, E c = 150 keV. 
Thermal case is plotted by a solid line, the non- 
thermal beam fluxes are £ - 10 s (dot), 10 9 (dash), 
10 10 (dash-dot), and 10" ergcirr 2 s _I (dash-dots). 



a result, there is only a small number of low- 
energy protons even for the high initial beam 
fluxes. 

The electron (proton) densities of the chro- 
mosphere obtained as a result of the NLTE cal- 
culations are of the order of 10 12 crrT 3 in the 
regions of the Ha centre and near wings for- 
mation (Figure H}. Thes e densities were com- 
puted using the code of iKasparova & Heinzell 
(2002) adapted f or the proton beams stu dies. 
It was shown by Bommi er et al.l (Il986bl) that 
depolarization of the hydrogen levels by col- 
lisions with the background electrons and pro- 
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Fig. 5. Hey line intensity profiles for the same set of 
models as in the Figure[4] 
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Fig. 6. Emergent fractional linear polarization Q/I 
profiles (%) computed close to the limb (fi = 0.11) 
for the same set of models as in the Figure|4] A pos- 
itive sign of Q/I means the tangential direction of 
polarization, while a negative sign is the radial one. 



tons becomes important even at the densities of 
the order of 10 10 cirT 3 . It is therefore expected 
that the depolarization effect will play a crucial 
role in the chromospheric conditions. 

The theoretical line profiles calculated for 
a set of the proton beam energy fluxes are plot- 
ted in the Figures [3] and [6] As it can be seen in 
Figure|6] a polarization degree of the emergent 
Ha line close to the limb is extremely small. 
Moreover, its orientation is tangential and not 
radial as it was expected for the low-energy 
proton impacts. This polarization is mainly due 
to resonance scattering of the radiation, while 
the impact polarization effect is not seen. The 
degree of 0.02% is well below any measurable 
value to day. 



5. Effect of electronic return currents 

It seems that linear polarization of the Ha line 
due to slow proton beams is not a good di- 
agnostics tool even if very restrictive condi- 
tions on the beam anisotropy are postulated. 
Another phenomena which could explain the 
observed radial polarization are the return cur- 
rents (RC) associated with the electron beams. 
The detailed studies of the collective plasma 
processes show that the return currents which 
neutralize a huge electric current of the beam 
can be formed if several plasma condi tions in 
the at mosphere are fulfilled (Nor man & Smith! 
119781) . The atomic excitations caused by the re- 
turn current were shown to be more important 
than the effect of the b eam itself in some depths 
dKarlickv et aDl2004. 

iKarlickv & Henouxl d2002l) propose 
the anisotropical impacts of the slow (few 
deca-eV) RC as an explanation for the ob- 
served linear polarization of the Ha line. 
However, all these calculations neglected the 
effects of depolarizing collisions and polarized 
radiative transfer. We have studied a simple 
model of monoenergetic electron beam with 
the initial energy flux of 1 .2x 10 12 erg crrr 2 s 
composed of the electrons with the energy 
10 keV penetrating the Fl atmosphere. The 
electron-hydrogen cross-sections for the 
dipolar transitions have been calculated by the 
semiclassical me thod with momentum transfer 
(BomMei] 120061) . These calculations show 
that return current is locally able to produce 
a significant atomic polarization, but it is 
decreased by the depolarizing collisions and 
high intensity of radiation. If all these effects 
are taken into account together with radiation 
transfer, the emergent linear polarization 
degree in the centre of Ha is only about 0.25% 
in the radial direction. That is still one order of 
magnitude below the measured values. 



6. Conclusions 

A net impact polarization of an electron beam 
+ RC seems to be a more promising explana- 
tion of the observed linear polarization than the 
proton beam impacts. However, collisional de- 
polarization by background electrons and pro- 
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tons together with the increased radiation in- 
tensity can destroy most of the impact atomic 
polarization. 

Several things have been simplified to 
make our models traceable: 

1. The semiempirical model Fl overestimates 
the atmosphereric temperature. As a result, 
depolarization and lines intensities are also 
overestimated. 

2. No time dependence was considered. 
Detailed time-dependent modeling could 
possibly explain some aspects of polariza- 
tion generation. 

3. The plane-parallel model has a limited ap- 
plicability and more complex geometry 
could be considered to interpret the obser- 
vations. 

4. If a RC is carried by a substantial num- 
ber of background electrons, they do not 
contribute to the collisional depolarization 
in the same way as the thermal electrons. 
That was not taken fully into account in our 
models. 

5. The beams are never purely directional and 
vertical and their scattering may lead to 
further decrease of the impact polarization. 
The monoenergetic electron beam + RC is 
not a satisfactory physical model. 

A quantitative estimation of all these effects 
should be a subject of the future studies. 
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